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DECLARATION UNDER 37C.F,R, 1.131 

We, Lav IVANOVIC, Paul FILSETH, and Mario GARZA, hereby declare that: 

1 . We are inventors of the subject matter recited in the claims of the above-identified 
application. 



2. Prior to October 9, 2001, we conceived of the idea of an automatic method to calibrate a 
masking process simulator that analyzes the difference between the aerial image produced by 
a simulator and the actual pattern produced by the masking process in order to improve the 
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results of the simulator, as described and claimed in our ideation. 

3. We conceived of the invention and reduced it to practice wbile walking for LSI Logic 
Coipoialion, having aprincipal place of business in Milpitas, California, in the OPC group, 
xnWch inchides LSI Logic peraonnel fiom CaUfiJinia as weU as from LSI Logic. International 
located in Russia. 

4. Attached Exhibit A is a copy of a documoit txanamittedby fiwsimile from LSI lo^o 
International to LSI logic (CA) on January 30, 2001, describing the •'Montenegro project 
(see pg. 3). Page 4 is a flow diagram of the project, and Hie blocks shown therein aie 
described on page 5-8 (see, for example, biockB12 and accompanying text). 

5. Attached Exhibit B is a screenshot of a directory structure stored on LSI logic computers 
showing "C* files in pr©-build state that were used implement the present inventiqa. As 
shown, all the files have dates of completion prior to February 29, 2000. 

6. Attached Exhibit C is a source code listing of the relevant portions of the source files for 
implementing the claimed operations above. 

I herd>y declare that all statements made herdn of my own knowledge are true and that 
all statements made on information and belief are believed to be true; and frirdier that these 
statemeats were made with flie knowledge that willful ftlse statements and the like so made are 
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pmiishaWe by fine or impiiswmuait, or both, wder Section 1001 of Title 18 of the United States 
Code and that such wiliafiil &lse statements may jeopardize the validity of the application or any 
patent issued fiioeon. 
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IlpoeKT «MoHTCHerpo» 
C.B.AjieniHH Oner 3a «HBapb2001 



3a npoiucAuiMii nepHoa rpynna npoanajmsHpOBaJia 
cocTOJiHHe OPC-iipoeKTa,HaMeTHJia nyxH pememui h nojiyHHna 

t KOHKpeTHUe peayJlbiaibl B HCKOTOpiWX HSUpaBWHHXX. 

BbiJio onpezie;ieHO,HTO flajiLHeHiuee npoABHSKCHHe CBxaano c 
pemeHueM saaai b Tpcx rjiaBHWX npo6jieMax, 

1 . Heo6xoAMMO flajibHefimee pasBHXHe h coBepmeHCTBOBaHHe 
nporpaMMM «MoJio-roB», to ecTb nporpawMW pcAaKrapOBaHiw 
flH3ajBMOB H CHMyjwuKH MacoKt PxA HCflOCTaxKOB flcflcTByiomefl 
BepcHH BHaBHJica npa oCpadcrrKc 6ojii.iiihx AH3aHHOB,Bce ohh 
coo6meHbi TTojiy aim HcnpaBJieHM«.Tpe6oBaHMa k «MojioTOBy» 
BoapacxaioT h b qbxsja c tcm , sxa nporpaMMa AOJi>KHa cxaxb 
HacTbK) eAHHoro nporpaMMHoro KOMiuieKca OPC. 

2. Heo6xOAHMO AajiiHeHiuee pa3BMTHe nporpaMMM OPC ,hto6w 
o6ecneHHTi> 6ojiee 6fcicTpyio o6pa6oTKy 6ojii»niHx AHsauHOB h, 
c Apyroii CTOpOHLi, eme OoJiee BMcuKoe KanecTBo o6pa6oTKH 
Majiwx nJiOTHBix AHsafiHOB (naMXTu). OjmoBpeMeHHo Heo6xoAHMO 
BecTH paapaOoncy nonbsoBaTejibCKoro wHTcp^jeHca. 

3. Heo6xoAHMO npHCTyiiHTB K paccMOxpeHMK) npo6jieMM aBTOMaxHsa- 
UHH npouecca KajiH6pOBKH cMwyjiaxopaB 3aBHCHMOCT!a ox cbohctb 
pe3Hcxa M onxHHecKott cHcxeMU. 

Mu pacHMTbiBaeMjHXO IIoji o6ecneHHx peajnoamuo n.l, 
flacx HaM onHcaHHe n,3 (KajiM6poBKM) he AocraxoHHo bmcokom 
ypOBHe (j|)opMajiH3auHH,a xaiOKC iqjHMex y^acTHe b o6cy3KAeHMM 
H peajiMsauHH noJibSOsaxejibCKoro uHxept^eiica. 

B luiawe peaimaaiwui n.2 6biJiH nonyneiai cjieAyiomHe 
pesyjibxaxw. 

E.EropoB rjiaBHoe BHHMaHHe yAeJiHJi 6jiOKy bwxoahoh 
MH(J)opMa^MH. CxaBHJiaci, aaAana nojiynirrb MaKCHManbHO 
KownaKTHoe onMcaHHC peayJibXHpyiomeH Macioi 6ojibmoro 
AHsaHHa, Pflfi 3xoro paspaGoxan cneuaaJibHUH $0pMaT,Ha3BaHHWii 



EXHIBIT A 
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OPC-4)opMaioM,Koropw*i iioaBOJiHci jiojiynaTb onMcaHHJi,o6iieM 

KOIOPWX Ha nopKflOK MCHbUJe TCK KOlOpwe MCnOJlbSOBaJlHCb flO CHX 

nop. 

yBeJlHHCHHe CKOpOCTH o6pa60TKH 6OJII1IUHX AHSaitHOB 

B03M0>KH0 Ha nyxH 6ojiee aKXHBHoro HcnojiL30BaHiM 
cMHTaKCMHecKMX npoucflyp H BUfleneHHH b j3H3attHe cxpyKxyp 
ciieuMajibHoro ama - npMMepoM tukhx cxpyKxyp Moryx 6i,iti, 
HEOTHCH, 6yKBLi, nepHOflHHecKHe cTpyKTypw,cTpyKTypw nonycKaio- 
mwe oflHOMepHyio o6pa6oTKy, onxHHecKH HeaaBHCHMHe noACTpyKiypbi 

Oahh h3 nepBWX (HanaJifcHMx) pesyjibxaxoB b 3xom HanpaBne- 
HMH nojiyHHJT C.PoflHH - iipoue^iypa BKwejieHiw hhcck xnna «6yKBa». 

M.MeflBefleBa peajiMSoeajia npoueflypy,no3Bon5iiomyK) noJiyHaxi, 
rpaHHuy H3o6paMceHHH hc coBnajiaioii^yio c rpaHHnefit HcxoAHoro 

JlH3aHHa,a C HeKOXOpbIM,3aAaHHBIM C^JBHrOM.STO flaeX BO3M0)KHOCXb 

peajiH30BaTb H^eio ^eHceHa o «paciiiHpeHHH» Hcxoj5Horo flHsaftwa. 
r.BejioKonbixoB npoBCJi cepmo SKcnepHMeHXOB no bmhwcjic- 

HHK) AByXCJlOflHblX ((taSOBblX MaCOKjBHHHCJWeMMX HaOCHOBC 

BsaHMOflelicxBHH flByx cjioeB nana (Poii n Island). I^ach TaKHX 
MacoKdwJia npeflno«eHa ;[^eHceHOM. OKaaajiocL, hto xaxHe mrckh 
^leilcTBHTejiBHo noBHmaioT KpyrHsny npo(|»wiH HHTeHCMBHocxH b 
pawoHe rpaHHUBi ,a cjie^OBaxejibHO Moryx yjiyHniHXb peayjitxax 
(|)oxojTHxorpa(J)HH.B to »e BpeMH B nJioxHE.ix sHsaiiHax co 
cjiojKHbiMH K0H4)HrypauHH Moryx BOSHHKaxb npo6jieMbi HexaaxKH 
CBo6oiiHoro npocxpaHCTBa. 

H.Ahhchmob npeflJiOMCHn o6myK> cxcMy HHxep4)eHca 
nojibsoBaxejw OPC-KOMmieKca.3xa cxeMa HanpaBnena Tlony wiR 
o6cy)KAeHHH. 

A.HepHyuiHH npHciynoJi k HsyneHHio KO^a nporpaMMW OPC. 
Uepen HHM nocxaBJicHa aaAana npeufloyianh nonpaBKn,KOTopMe 
yjiyMiiiai xecxHpyeMOCTb nporpaMMW a xaicace JiHKBHAHpyiOT 
HMeiomHecH ceftnac 3(t$eKXbi «yxeHKH naMaTH». 
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LSI Logic International 

MONTENEGRO project MecaiHHbia oT Hex 

HcnoJiHHTejib; H.^.Ahhchmob IlepHoj.: gnaapb 2001 r. 

OcHOBHue pesyjibxaTU. 

1 Pa3pa6oTKa ajiropHTMOB. 

• pa3pa6oTaHa o6maji SjioK-cxcMa. MHJt«p4>«8ca noJibsoBaTejisi ajih pc- 
uieHHfl 3aAaMM KoppeKU»iH; 

• onwcwo <|)yHKBtHOHHpoBaHMe 0Tjae-ii»in»ix 6jiokob AaHHott cxeMu; 
BjioK-cxewa h onHcaHMe 4)yHKiwoHHpoBaHHa ee 6/iokod npHJiaraexcH. 

2 IIjiaHW Ha cjie/iyioiii^iii MecHU. 
2.1 Pa3pa6oTKa aJiropHTMOB. 

• npoflOJiJKCJtiHe pa3pa6oTKH HMTcp+ettca nojib30BaTe^>i; 

• nporpawMHpoBaHHe oTOejibHwx 6jiokob HHTep<i)elica; 
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1 Flowchart of user interface and 

description of its links and function- 
General flowchart of user interface is represented at fig-L 
Description of its action is given below. 

1.1 Block Bl: composition of task and setup of 
general parameters. 

• determination of paths of search for 

- input data; 

- parameter files; 

- output data. 

• determination of general parameters of design: 

- input design gds-file name; 

- top structure name in input gds-file; 

- format of output data; 

- output file name; 

- description of layers for input gds-file; 

- description of layers for output mask file; 

- types of mask used. 

• choice of models and their parameters setup: 

- optical model; 

- simulator model; 

- resist and exposure model; 

- etching model. 

• determination of the type of design problem; 

- simple design; 

- complex design. 



1 
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1.2 Block B2 : setup of big design parameters. 

Setup of the specific parameters Tor processing of big design, such as: 

• classification algorithm parameters; 

• optimization algorithm parameters; 

• mask composition algorithm parametersj 

• parameters of quality check and additional correction. 

1.3 Block B3 : setup of periodic design parameters. 

Setup of the specific parameters for processing of periodic design, such as: 

• classification algorithm parameters; 

• optimization algorithm parameters; 

• mask composition algorithm parameters; 

• parameters of quality check and additional correction. 

1.4 Block B4: setup of small design parameters. 

Setup of the specific parameters for processing of small design, such as: 

• classification algorithm parameters; 

• optimization algorithm parameters; 

• mask composition algorithm parameters; 

• parameters of quality check and additional correction. 



2 
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1.5 Block B5: design transformer. 

Completes decomposition of initial complex design on to parti^tl fragments, 
which have move simple structure (periodic memory, ^^^^^'/^^^^^f^f ' ^^^^^^ 
Generation of the data flow for further processing of fragments and formation 

of the task stack. 

Design fragmentation may be completed: 

• by cell type or list of cell types-, 

• by cell name or list of cell names: 

. by chosen region of list of regions (or with exclusion of some regions); 

• by density chart. 

1 6 Block B6: classification, optimization, local mask 
composition, quality check up and additional cor- 
rection for big design. 

This block contains basic OPC procedures modified for «>inplex designs The 
detLd description of these procedures will be added m nearest future. 

1,7 Block B7: statistics. 

This block will contain a set of procedures, which enable to user: 

• to obtain numerical and qualitative estimation of the design (volume, 
complexity etc); 

• to display and analyze the sequense of stages of design processing; 

• to create protocols of design correction; 

. to search the information current status of processes and on resources 

in use; 

• to create different log-files. 



3 
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1.8 Block B8: visualization. 

The rocedures of this block enable user to execute such operation as: 

• display information given by block B7 " Statistica" ; 

• to rnn visualization tools, such as Molotof or Artwork Editor, for scan 
and edition of input designs, masks, simulation results, etc; 

• to build hystogxams and graph which illustrate results of the design 
processing; 

• to convert input data represented in different graphics formats(gds, bgl, 
tgl, MEBES, CIF, etc). 

1.9 Block B9: mask composition, quality check up 
and additional correction for periodic design. 

This block contains basic OPC procedures modified for periodic designs. The 
more detailed description of these procedures will be added in nearest future. 

1.10 Block BIO: mask composition, quality check up 
and additional correction for small design. 

This block contains basic OPC procedures modified for small designs. The 
more detailed description of these procedures will be added in nearest future. 

1.11 Block Bll: mask composer. 

This block executes the procedure of output mask composition for entire 
design from local and fragment masks and forms output gds-file. 

1.12 Block B12: calibration. 

This block will be intended to make (in automatic and/or interactive mode) 
the procedure of adjustment of models parameters to empirical data obtained 
for some test patterns. Such adjustment will give more adecuate description 
of phisical and technological processes (diffraction image formation, exposure, 
post-exposure development and mask making) for ope 
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Sxibclaims 1-c and 1-e semedge.c May 25, 1999 



This file contains a computer-vision algorithm for detecting the 
edges of features on a digitized gray-scale image of a pattern, 
and a procedure for drawing the detected edges onto the digital 
image . 



* 

* Find the features in a SEM 
★ 



ttinclude "Screens . h" 

#define ZOOM_STEPS 20 

#define HASH_SIZE 4096 

#define TOO_WIDE 100 // Limit on SEM feature edge width 

#define HIST_SIZE 20 // HIST_SIZE > 2 * log2 (TOO_WIDE) 

#define NEAR{v, t, r) ((unsigned) { (v) - (t - r) ) <= (r * 2)) #define 

PIF if (cfp) fprintf{cfp, 

#def ine ENDPIF ) ; 

ttdefine ENDPIFF ), fflush(cfp) ; 

typedef struct { 

char which_dir; 
} connectionsrec, *connections ; 

typedef struct locREC { 

struct locREC *next; 

int color; 

short X, y; 

short feature_id; 

connect ionsrec left, right; 

char cleft, cright; 

char reached; 
} locrec, *loc; 

typedef struct pairingREC { 

struct pairingREC *next; 

struct splotchREC *sl, *s2; 

int count; 
} pairingrec, *pairing; 

typedef struct splotchREC { 
struct splotchREC *next; 
struct splotchREC *mesa; 
int id; 
int dir; 
loc *edge; 
pairing partners; 
int polarity; 



EXHIBIT C 



} splotchrec, *splotch; 



^•k1fk'kifk1tit1t1t1t'k**'kiticirifk*'k*itit*i(1t Variables 
***«********4r*4r*ilr*** ************ y 

Static int bright [256] ; 

static int variation; 

static int lo_x, hi_x, lo_y, hi_y; 

static int curcolor = 0; 

static int array_size = 0; 

static int *profile, *al, *a3, *a9, *rising, ♦fallings- 
static Display *theDisp; 
// static FILE *cfp = (FILE *) -1; 
static FILE *cfp = NULL; 
static loc lochash [HASH_SIZE] ; 
static int grad_radius, edge_width, box_size; 
static int features [500] ; 
static int tiistogram [HIST_SIZE + 1]; 
static int feature id; 



static connect ions rec empty_connections = {-l}; 
extern Colormap my_colormap; 

y******************* ********* Functions 
****************** ***************^ 

#define EOFF 3 

draw_feature_edges_on_SEM (Widget w, grapliics_data *data) 

{ 

splotcll S; 
loc p; 

int i, ix, iy, j, ox, oy, bad, obad, phase; 
Display *dpy = XtDisplay (w) ; 

Window win = XtWindow (w) ; 

GC gc; 

pmversion pm; 
double factor; 
connections c, r; 

static string color [] = {"Violet", "Green", "Red", "Blue", "Yellow"}; 



pm = data->curpm; 

factor = {(double) pm->zoom_level) / ZOOM_STEPS; 

gc - data->gc; 

gc = data->gcSEMLabel; 

for (j = 0; j < 5; ++ j ) { 

XSet Foreground (dpy, gc, color_code (color [ j ] , w) ) ; 

loop_over_list (s, data->f eatures) { 
obad = ox = oy = -1; 

for (i = 0/ s->edge[i] != NULL; ++i) { 
p = s->edge [i] ; 

ix = (int) (factor * p->x + 0.5); 
iy = (int) (factor * p->y + 0.5); 
if (data->flipX) 

ix = data->sem->xres - ix; 



} 



} 



} 



} 



if (data->f lipY) 

iy = data->sem->yres - iy; 
/* Alternate color rules 
phase = s->mesa->id % 5; 

if (s->polarity >= 2) 

phase = 2; 
else if (s->polarity <= -2) 

phase = 3 ; 
else if {s->polarity == 0) 

phase = 4; 

else 

phase = 1; 

*/ 

if {s->polarity >= 2) 

phase = 4; 
else if (s->polarity <= -2) 

phase = 3; 

else 

phase = -1; 
if (phase == j) { 

// XDrawPoint (dpy, win, gc, ix + EOFF, iy + EOFF) ; 
if (i 1= 0) 

XDrawLine (dpy, win, gc, 

ix + EOFF, iy + EOFF, ox + EOFF, oy + EOFF) ; 

} 

ox = ix; 
oy = iy; 
obad = bad; 



static int brightness (graphics_data *data, int x, int y) 

{ 

pmversion pm; 
windo sem; 
int c, pix; 
XColor color; 

pm s= data->curpm; 
sem = data->sem; 
if (x < 0) 

X = 0; 
if (x >= sem->xres) 

X = sem->xres; 
if (y < 0) 

y = 0; 
if (y >= sem->yres) 

y = sem->yres; 
c = pm->explmage->data [y * pm->width + x] & 255; 
if (bright [c] == 0) { 

pix = (int) XGetPixel (pm->explmage, x, y) ; 

color .pixel=pix; 

XQueryColor (theDisp, my_colormap, &color) ; 
bright [c] = (int) 




((color. blue + color. green + color. red) * 0.333333 / 

65.536) ; 

if (bright [c] == 0) 
bright [c] = 1; 

} 

return (bright [c] ) ; 

} 

static expand_arrays (sem) 
windo sem; 

{ 

int needed; 

needed = (sem->xres > sem->yres ? sem->xres : sem->yres) + 25; 
if (needed <= array_size) 

return; 
array_size = needed; 
al = newarray (int, needed) + 11; 
profile = al; 

a3 = newarray (int , needed) + 11; 
a9 = newarray (int, needed) + 11; 
rising = newarray (int, needed) + 11; 
falling = newarray (int , needed) + 11; 

} 

static prof ile_line {graphics_data *data, int start, int dx, int dy) { 
int X, y, pix, i, j; XColor color; int val, last, diff; Boolean 
falling; windo sem; 

sem = data->sem; 
expand_arrays (sem) ; 
if (dx 0) { 
y = 0; 
X = start; 
} else if (dy == 0) { 
X = 0; 
y = start; 
} else if (dx == 1 && dy == 1) { 
if (start >= 0) { 
X = 0; 
y = start; 
} else { 

X = - start; 
y = 0; 

} 

} else if (dx == 1 && dy == -1) { 
if (start >= 0) { 
X = 0; 

y = sem->yres - 1 - start; 
} else { 

X = - start; 

y = sem->yres - 1; 

} 

} else { 
return; 

} 

last = 0; 



falling = True; 

for (i = 0; True; x += dx, y += dy, ++i) { 
if (x < 0 I I X >= sem->xres) 
break; 

if (y < 0 I I y >= sem->yres) 
break; 

val = brightness (data, x, y) ; 
profile [i] = val; 

} 

profile [i] = -1; 

} 

static int inf o_content {graphics_data *data, 

int start, int dx, int dy, Boolean calibrate) 

{ 

int i, last, v, a, medium, total; 
int low, high; 
double fl, f2, f3; 

prof ile_line (data, start, dx, dy) ; 
medium = total = 0; 
low = high = 0; 

for (i = 3; prof ile [i] >= 0; ++i) { 
last = profile [i - 3]; 
V = profile [i] ; 
a = V - last; 
if (a < 0) 

a = - a; 
if (a > 80) 

a = 0; 
total += a; 
if (v < 200) 

++ low; 
if (v > 800) 

++ high; 

} 

if (calibrate) { 

variation += total / i; 
return (0) ; 

} 

fl = (double) high / i; 
f2 = (double) low / i; 
f3 = ( (double) total / i) / variation; 
if (fl > 0.9 && f3 < 0.5) 
i = 1; 

else if (f2 > 0.4 && f3 < 0.8) 

i = 2; 

else if (fl + f2 > 0.9 f3 < 0.6) 

i = 3; 
else 

i = 0; 
return ( i ) ; 

} 

static int f ind_border (graphics_data *data, int dx, int dir) 

{ 

windo sem; 



int out, in, mid, dy, range, last; 



sem = data->sem; 
dy = 1 - dx; 

range = dx ? sem->yres : sem->xres; 
if (dir == 1) { 

out = -1; 

in = range / 2; 
} else { 

out = range; 

in = range / 2; 

} 

last = -1; 

while (in != out + dir) { 
mid = (out + in) / 2; 
if (mid == last) 

printf { "Screwup\n" ) , exitd); 
if (info_content (data, mid, dx, dy. False) == 0) 
in = mid; 

else 

out = mid; 
last = mid; 

} 

return (in) ; 

} 

static identify_border_areas (graphics_data *data) 

{ 

windo sem; 



sem = data->sem; 
variation = 0; 

info_content (data, (int) (sem->xres * 0.4) 

inf o_content (data, (int) (sem->xres * 0.5) 

info_content (data, (int) (sem->xres * 0.6) 
variation /= 3; 

lo_y = f ind_border (data, 1, 1); 

hi_y = f ind_border (data, 1, -1) ; 

lo_x = f ind_border (data, 0, 1); 

hi_x = f ind_border (data, 0, -1); 



0, 
0, 
0, 



1, 
1, 



True) ; 
True) ; 
True) ; 



PIF "Border: 
ENDPIFF } 



X = %d-%d, y = %d-%d\n\n", lo_x, hi_x, lo_y, hi_y 



static void record_j)air_stats (int width) 

{ 

int ix, pow; 



if (width <= 3) { 

ix = width; 
} else { 

for (pow = 2; (1 << pow) <= width; ++pow) ; 
ix = (width >> (pow - 2)) + pow * 2 - 4; 

} 

PIF "Width = %d, ix = %d\n", width, ix ENDPIF 
++ histogram [ix] ; 



static int ix_to_n(int i) 

{ 

return ({(2 + (i & 1)) << (i » l) ) >> l) ; 

} 

Static void analyze^pair^statistics () 
int i, V, max, ix; 
ix = max = 0; 

for (i = 0; i < HIST_SIZE; ++i) { 
if (histogram [i] != 0) 

PIF "%d: range = %d to %d, num of pairs = %d\n", i, 
ix_to_n(i), ix_to_n(i +1) - 1, 
histogram [i] ENDPIF 
V = (histogram[i] + histogram[i + 1] ) * i; 
if (v >= max) { 
max = v; 
ix = i; 

} 

} 

edge_width = ix_to_n(ix + 1) ; 

PIF "Typical width of suspected edges = %d pixels\n", edge_width 
ENDPIF 

printf ("Typical width of suspected edges = %d pixels\n", 
edge__width) ; 

grad_radius = edge_width / 4; 
if (grad_radius < 1) 

grad_radius = 1; 
if {grad_radius > 5) 

grad_radius = 5 ; 
// box^size = edge_width * 5 / 2/ 
box_size = edge_width * 2; 

} 



static int f ind_f eatures (int lo, int hi) 

{ 

int i, e, last, base, j, min, count; 

j = (al[lo] + al[lo + 1]) » 1; 

for (i = lo - 10; i < lo; ++i) 

al[i] = j; 

j = (althi] + al[hi -1]) >> 1; 



for (i 


= hi + 10; i 


> hi; --i) 




al[i] 


= j; 






for (i 


= lo - 9; i < 


= hi + 9; ++i) 




a3[i] 


= (al[i - 1] 


+ al [i] + al [i 


+ 1]) 


a3 [lo - 


10] = a3 [lo 


- 9]; 




a3 [hi + 


10] = a3 [hi 


+ 9]; 




for (i 


= lo - 7; i < 


= hi + 7; ++i) 




a9[i] 


= (a3[i - 3] 


+ a3 [i] + a3 [i 


+ 3]) 



for (i = lo - 2; i <= hi + 2; ++i) { 

rising [i] = profile [i] - prof ile [i - 1] + a3 [i + 1] - a3 [i - 1] + 



a3 [i + 1] - a9 [i - 5] ; 
falling [i] prof ile [i] - prof ile [i + 1] + a3 [i - 1] - a3 [i + 1] 

a3 ti - 1] - a9[i + 5] ; 

} 

count = 0; 
last = -99; 

for (i = lo/ i <= hi; ++i) { 
e = rising [i] ; 

if (e > 200 && e > rising [i - 1] && e >= rising [i + 1] && 
e > rising [i - 2] && e >= rising [i + 2]) { 
// fprintf(cfp, " r[%d] = %d\n", i, e) ; 
last = i; 

} 

e = falling [i] ; 

if (e > 200 && e > falling [i - 1] && e >= falling [i + 1] && 
e > falling [i - 2] && e >= falling [i +2]) { 
// fprintf(cfp, " f[%d] = %d\n", i, -e) ; 
if (last != -99) { 
base = a3 [last - 2]; 
if (base < a3 [i + 2] ) 
base = a3 [i + 2] ; 
if (i <= last + 1) { 

min = prof ile [prof ile [i] < prof ile [last] ? i : last]; 
} else { 

min = a3 [last + 1] ; 
for (j = last + 2; j < i; ++ j ) 
if (a3 [j] < min) 
min = a3 [j] ; 

} 

PIF "Pair at %d-%d, edges = %d, height = %d\n", 
last, i, 

rising [last] + falling [i], 
min - base ENDPIF 
if (i - last > TOO_WIDE) { 

PIF " Rejected: too wide\n" ENDPIF 
} else if (min - base <= 0) { 

PIF " Rejected: too low\n" ENDPIF 
} else { 

record_j5air_stats (i - last); 

features [count++] = last; 

features [count++] = i; 

features [count] = -99; 

} 

last = -99; 

} 

} 

} 

if (cfp) fflush(cfp); 
return (count) ; 



static int *cross_section (graphics_data *data, int start, int dx) { int 
dy; int i, x, y, n, dir; loc p; int *ptr; 



dy = 1 - dx; 

if (cfp) fprintf(cfp, "cross_section %d, dx=%d, dy=%d\n", start, 
dx, dy) ; 

prof ile__line (data, start, dx, dy) ; 
if (dx == 1) 

n = f ind_features (lo_x, hi_x) ; 
else 

n = f ind_features (lo_y, hi__y) ; 
ptr = newarray (int , n * 2 + 1) ; 
for (i = 0; i < n; ++i) { 
if (dx == 1) { 

X = features [i] ; 
y = start; 
} else { 

X = start; 

y = features [i] ; 

} 

ptr (2 * i] = x; 
ptr [2 * i + 1] = y; 

} 

ptr [2 * n] = -99; 
return (ptr) ; 

} 

static int octant, mag, level; 
static int delta_x, delta__y; 
static string dirname[] = { 

"ENE", "NNE", "NNW", "WNW" , "WSW", "SSW", "SSE", "ESE" }; 

/* This approximately converts cartesian coordinates to polar 
coordinates, 

based on pretending an octagon is a circle. It returns the angle in 
45-degree chunks and leaves the radius in "mag". */ 

static int cartes ian_to_polar (dx, dy) 

int dx, dy; 

{ 

static int octant_map [8] - {O, 1, 3, 2, 7, 6, 4, 5}; 
int adx, ady; 
int ix; 



if (dx < 0) 

ix = 2, adx = - dx; 
else 

ix = 0, adx = dx; 
if (dy < 0) 

ix += 4, ady = - dy; 
else 

ady = dy; 
if (adx < ady) 

++ ix, adx >>= 1; 
else 

ady >>= 1; 
mag = adx + ady; 
return (octant_map [ix] ) ; 



static splotch new_splotch (splotch *plis) 
{ 

splotch s; 

s = new(splotchrec) ; 
s->id = feature_id; 
s->next = *plis; 
*plis = s; 
s->mesa = NULL; 
s->partners = NULL; 
s->edge = NULL; 
s->dir = 0; 
s->polarity = 0; 
return (s) ; 

} 

#define EDGE_EFFECT 5 

static f ind_gradient (graphics_data *data, int x, int y) 
{ 

static int octant_map [8] = {O, 1, 3, 2, 1, 6, 4, 5}; 
int i, j, k, dx, dy, adx, ady; 
int val [9] ; 

for (k = 0, i = -grad_radius; i <= grad__radius ; i += grad_radius) 
for (j = -grad_radius; j <= grad_radius; j += grad_radius) 
val[k++] = brightness (data, x + i, y + j); 
dx = (val [6] - val[0] + val [7] - val [1] + val [8] - val [2]) / 3; 
dy = (val [2] - val[0] + val [5] - val [3] + val [8] - val [6] ) / 3; 
if (x <= lo_x + EDGE_EFFECT) { 
if (dx < 0) 
dx = 0; 

} else if (x >= hi_x - EDGE_EFFECT) { 
if (dx > 0) 
dx = 0; 

} 

if (y <= lo_y + EDGE_EFFECT) { 
if (dy < 0) 
dy = 0; 

} else if (y >= hi_y - EDGE_EFFECT) { 
if (dy > 0) 
dy = 0; 

} 

dy = - dy; 
delta_x = dx; 
delta_y = dy; 
level = val [4] ; 

octant = cartesian_to_polar (delta_x, delta_y) ; 



static init_hash_table 0 

{ 

int i; 

for (i = 0; i < HASH_SIZE; lochash[i++] = NULL); 

} 



static int hcalls, hsteps, hmax, hpoints, hrepeats; 



static loc lookuploc(x, y) 
int X, y; 

{ 

int ix, count; 
loc p, *ptr; 

++ hcalls; 

ix = (x + y) + (x << 3) - (x >> 5) + (x << 6) - (x >> 2) ; 
ix &= HASH_SIZE - 1; 
ptr = lochash + ix; 
count = 0; 

while (*ptr != NULL) { 
p = *ptr; 

if (p->x == X Sc& p->y == y) 

return (p) ; 
ptr = &p->next; 
++ count; 

} 

if (count > hmax) 

hmax = count; 
hsteps += count; 
++ hpoints; 
p = newdocrec) ; 
*ptr = p; 
p->next = NULL; 
p->color = 0; 
p->X = x; 

p->y = y; 

p->f eature_id = -1; 
p->left = empty_connections; 
p->right = empty_connections; 
p->cleft - p->cright = 0; 
p->reached = 0; 
return (p) ; 

} 

static loc stack [5000]; 
static int sp; 

static loc outline [5000] ; 
static int stepdir [5000] ; 

static int step_x[8] = {l, 0, -1, -1, -1, 0, 1, l}; 
static int step_y[8] = {l, 1, 1, 0, -1, -1, -1, O}; 
static int absangle[8] = {O, 1, 2, 2, 4, 3, 2, l}; 

#define B 4 
int pgf = 0; 

static loc follow_gradient ( 

graphics_data *data, int x, int y, Boolean to left, int cut) 

{ 

int origdir, recentdir, angle; 

int oc, i, j, k, newx, newy, curx, cury, max, maxoc, delta, steepness 
int lastoctant, dir, x_change, y_change, rightangle, prevmax; loc p, 
r; short *ptr; Boolean looped; splotch s; int levels [4] , level_sum; i 



delx[B + 1], dely[B +1]; connections c; char *connec; Boolean foo; 
Boolean boundary; 



if {to_left != 0 6c& to_left != 1) { 
print f { "Screwup\n" ) ; 
exit (1) ; 

} 

++ feature_id; 

PIF "%d: Tracing %c from %d, %d\n'\ 

f eature_id, "RL" [to_left] , x, y ENDPIFF 

curx = x; 
cury = y; 
looped = False; 
k = 0; 

level_sutn = 0; 

for (i = 0; i <= B; ++i) 

delx[i] = dely[i] = 0; 
for (i = 0; i < 4; ++i) 

levels [i] 0; 

p = lookuploc (curx, cury); 
outline [0] = p; 

rightangle = to_left ? 2 : -2; 

c = to_left ? &p->left : &p->right; 

lastoctant = -1; 

boundary = False; 

prevmax = 50; 

for (j = 0; True; ++ j ) { 

++ pgf; 

f ind_gradient (data, curx, cury); 
steepness - mag / 10; 

level_sum += level - levels [j & 3]; 
levels [j & 3] = level; 

delx[B] += delta_x - delx[j & 3]; 
dely[B] += delta_y - dely[j & 3] ; 
delx[j & 3] = delta_x; 
dely[j & 3] = delta_y; 

recentdir = cartesian_to_polar (delx [B] , dely[B]); 
if (j >- 4) 

steepness = mag / 40; 

origdir = octant; 

if (c->which_dir < 0) { 

max = -1; 

maxoc = origdir; 

dir = (recentdir +2) & 7; 

for (i = 0; i < 4; ++i) { 
/* 

if (i 3 && j > 2) 
break; 

*/ 

if {to_left) 

oc = (origdir + i) & 7; 



else 

oc = (origdir - i - 1) & 7; 
newx = curx + step_x[oc]; 
newy = cury - step_y[oc]; 
if (j > 0) { 

q = outline [j - 1] ; 

if (newx == q->x && newy == q->y) 
continue; 

} 

f ind_gradient (data, newx, newy) ; 
if (j >= 3) { 

angle = absangle [ (octant - recentdir) & 7] ; 
if (angle > 1) { 
if (angle > 2) 

continue; 
mag >>= 1; 

} 

delta = (level_sum >> 2) - level; 
if (delta > 0) 

mag -= delta; // + (delta >> 1); 

} 

if (mag > max) { 
max = mag; 
maxoc = oc; 

} 

} 

} else { 

maxoc = c->which_dir; 

} 

stepdir[j] = maxoc; 

x_change = step_x [maxoc] ; 

y_change = - step_y [maxoc] ; 

curx = curx + x_c]iange; 

cury = cury + y_cliange; 

lastoctant = maxoc; 

p = loolcuploc (curx, cury); 

outline [j + 1] = p; 

c = to_left ? &p->left : &p->riglit; 

if (j < 15 I I j % 10 == 0) { 

if (cfp) fprintf(cfp, "%s%s %d (%d,%d)", 
(Jc % 4 3) ? "\n\t" : " 
dirname [ (origdir + riglitangle) & 7], 
steepness, curx, cury); 

++ k; 

) 

if (p->f eature_id == -i - feature_id) { 
if {! looped) { 

PIF "Detected a loop at %d,%d\n", curx, cury ENDPIF 
looped = True; 

} 

} else if (p->f eature_id < 0) { 
if (j >= cut) 

p->f eature_id = -1 - feature_id; 
} else { 

if (c->wliich_dir >= 0) { 

// if (p->reached & (1 << to_lef t) ) { 

if (cfp) fprintf (cfp, "\nEntered previously taJcen patli\ 



if (j < 3) 

return (NULL) ; 
break; 

} 

++ hrepeats; 

} 

if (looped) { 

for (1=1; i <= j; ++i) { 
q = outline [ j + 1 - i] ; 
if (q == p) 
break; 

} 

if (cfp) fprintf (cfp, "\nEdge forms a loop of %d points\ 

if (i < 30 && j < 50) 

return (NULL) ; 
break; 

} 

if (max < 40 /* Sc& prevmax < 40 */) { 

PIF "\nEdge faded out (grad = %d) at %d,%d\n", 

max, curx, cury ENDPIF 

if (j < 50) 

return (NULL) ; 
break; 

} 

if (x_change <= 0 && curx <= lo_x + 2 | | 

x_change >= 0 && curx >= hi_x - 2 | | 
y_change <= 0 && cury <= lo_y + 2 | | 
y_change >= 0 && cury >= hi_y - 2) { 
if (j > 0) { 

if (cfp) fprintf (cfp, "\nEdge reached boundary\n") ; 
boundary = True; 
if (j < 2) 

return (NULL) ; 
break; 

} 

prevmax = max; 

} 

if (j > 14) { 

q = outline [ j - 13] ; 

if {(unsigned) (p->x - q->x +4) <= 8 && 

(unsigned) (p->y - q->y + 4) <= 8) { 
PIF "\nEdge not making progress\n" ENDPIF 
if (j < 50) 

return (NULL) ; 
break; 

} 

} 

if (j >= 4000) { 

if (cfp) fprintf (cfp, "\nJust keeps going and going... \n 
break; 

} 

} 

++ j; 

if (j < cut && ! boundary) 
return; 



if (cfp) fprintf(cfp, "Adding splotch %d\n", f eature_id) , 
fflush(cfp) ; 

s = new_splotch(&data->features) ; 
s->dir = to_left; 
s->edge = newarray ( loc , j + 2) ; 
if (j < 7) 

i = 0; 
else if (j < cut) 

i = j / 2; 
else 

i = cut ; 

for (k = 0; i <= j ; ++k, ++i) { 
p = outline [i] ; 
s->edge [k] = p; 

c = to_left ? &p->left : &p->right; 
if (i != j) { 

connec = to_left ? &p->cleft : &p->cright; 

*connec |= 1 << stepdir[i]; 

if (c->which_dir >= 0 c->which_dir != stepdir[i]) { 
if (cfp) fprintf (cfp, "Inconsistency at %d,%d: old=%d, 

new=%d\n" , 

p->x, p->y, c->which_dir, stepdir[i]); 
c->which_dir = stepdir[i]/ 

} 

if (k != 0) { 

connec = to_left ? &p->cright : &p->cleft; 
*connec |= 1 « {stepdir[i - 1] ^ 4) ; 

} 

if (i > 4 I I j < 7) { 

p->reached |= 1 << to left; 

} 

if {p->f eature_id < 0) 

p->f eature_id = feature id; 

} 

s->edge[k] = NULL; 

PIF "Stacking %d,%d\n", outline [j ]- >x, outline [j ] ->y ENDPIF 
stack [++sp] = outline [j]; 
for (i = 40; i < j; i += 20) { 
stack [++sp] = outline [i]; 

PIF "Stacking %d,%d\n", outline [i] ->x, outline [i] ->y ENDPIF 

PIF "\n" ENDPIFF 
return (outline [ j ] ) ; 

} 

static int bottom, top; 

static int f ind_range (char c, int go_left) 

{ 

int i, count; 

if (c == 0) { 

bottom = top = -1; 
return (-1) ; 

} 

count = 0 ; 



for (i = 0; i < 16; ++i) { 
if ((c >> (i & 7)) & 1) { 
if (count >= 4) 

break; 
count = 0; 
} else { 

++ count; 

} 

} 

if (i == 16) { 

bottom = top = -2; 
return (-2) ; 

} 

bottom = i & 7; 
top = bottom; 

for (i = bottom +1; i < bottom + 4; ++i) { 
if ( (c » (i & 7) ) & 1) 

top = i & 7; 

} 

return (go_left ? bottom : top) ; 

} 

Static int break_count; 

static sever (loc p, int goleft, int dir, int depth) 

int X, y, k; 
loc q; 

char *ptr, *qptr; 
Boolean foo; 

foo = cfp != NULL && p->x == 485 && p->y == 231; 
ptr = goleft ? &p->cleft : £cp->cright; 
if ((*ptr & (1 << dir)) == 0) 
return; 

if ((unsigned) (p->x - 493 +3) <= 6 && (unsigned) (p->y - 230 + 3) 
<= 6) 

if (cfp) fprintf (cfp, "Severing %d,%d: %c dir=%d, #=%d, 
depth=%d\n", 

p->x, p->y, "RL" [goleft] , dir, break_count, depth); 
*ptr &= ^ (1 << dir) ; 
X = p->x + step_x[dir]; 
y = P->y - step_y[dir]; 
q = lookuploc(x, y) ; 

if (foo) fprintf (cfp, "Cutting link to %d,%d (%x<->%x, 
reached=%d) \n" , 

q->x, q->y, q->cleft & 255, q->cright & 255, q->reached) ; 
qptr = goleft ? &q->cright : &q->cleft; 
*qptr &= ~ (1 << (dir ^ 4)); 
if (*qptr == 0 && q->reached != 3) { 
for (k = 0; k < 8; ++k) { 

sever(q, goleft, k, depth + 1); 

} 

} 

} 

#define LOOKAHEAD 10 



static int pick_best__path (loc p, int goleft, int preferred, int depth) 
{ int step, i, j, k, keeper, ix, x, y, going, choice, maxstep; int 
dirs[4]; int prevcolor [4 * LOOKAHEAD] ; loc path [4 * LOOKAHEAD] ; int 
route [4 * LOOKAHEAD]; int savecolor; loc q, r; char qcon, *ptr; Boolean 
foo; 

if (depth > 10) 
return ( -1) ; 

foo = NEAR(p->x, 111, 3) && NEAR(p->y, 321, 3); 
savecolor = curcolor; 

ptr = goleft ? &p->cleft : &p->cright; 
k = 0; 

if (goleft) { 

for (i =0; i < 4; ++i) { 

j = (preferred + i) & 7; 
if (*ptr fic (1 << j) ) 
dirs [k++] = j ; 

} 

} else { 

for (i = 0/ i < 4; ++i) { 

j = (preferred - i) & 7; 
if (*ptr & (1 << j) ) 
dirs [k++] = j ; 

} 

} 

going = k; 

PIF "pick_best_path %d,%d (%c) dirs=%d, depth=%d\n", 

p->x, p->y, "RL" [goleft] , going, depth ENDPIF 
for (i = 0; i < k; ++i) { 

ix = i; 

x = p->x + step_x [dirs [i] ] ; 
y = P->y - step_y [dirs [i] ] ; 
q = lookuploc(x, y) ; 
prevcolor [ix] = q->color; 
route [ix] = dirs [i] ; 
q->color = savecolor + ix; 
curcolor = q->color + 1; 
path[ix] = q; 

} 

for (step = 0; step < LOOKAHEAD - 1 && going > 1; ++step) { 
for (i = 0; i < k; ++i) { 

ix = i + (step +1) * 4; 
path[ix] = NULL; 
q = path [ix - 4] ; 
if (q NULL) 
continue; 

qcon = goleft ? q->cleft : q->cright; 
keeper = f ind_range (qcon, goleft); 
if (keeper < 0) { 
going; 
continue; 

} 

if (top != bottom) { 
if (q == p) 

return (-1) ; 

keeper = pick_best_j)ath(q, goleft, keeper, depth + 1); 



if (keeper == -1) 
return (-1) ; 

} 

X = q->x + step_x [keeper] ; 
y = q->y - step_y [keeper] ; 
r = lookuploc(x, y) ; 
if (r->color >= savecolor) { 

" going; 

continue; 

} 

prevcolor [ix] = r->color; 
route [ix] = keeper; 
r->color = savecolor + ix; 
curcolor = r->color + 1; 
patli[ix] = r; 

} 

} 

if (going == 0) 

step; 
maxstep = step; 

for (choice = 0; clioice < k; ++choice) { 
ix = choice + step * 4; 
if (path[ix] 1= NULL) 
break; 

} 

if (choice >= k) 

return (-1) ; 
*ptr = 1 << dirs [choice] ; 
curcolor = savecolor; 
for (i = 0; i < k; ++i) { 

for (step = 0; step <= maxstep; ++step) { 
ix = i + step * 4; 
q = path [ix] ; 
if (q != NULL) 

q->color = prevcolor [ix] ; 

} 

} 

for (i = 0; i < k; ++i) { 
if (i == choice) 
continue; 

for (step = 0; step <= maxstep; ++step) { 
ix = i + step * 4; 
q = path[ix] ; 
if (q == NULL) 

break; 
keeper = route [ix]; 

ptr = goleft ? &q->cright : &q->cleft; 

*ptr &= - (1 « (keeper ^ 4)); 

qcon = *ptr; 

// *ptr = 0; 

if (step != 0) { 

ix = i + (step - 1) * 4; 

q = path [ix] ; 

ptr = goleft ? &q->cleft : &q->cright; 

*ptr &= - (1 << keeper); 

// q->cleft = q->cright = 0; 

} 



if (qcon != 0) 
break; 

} 

} 

return (dirs [choice] ) ; 

} 

static int callno; 

static break_links (splotch s) 

{ 

int i, goleft, len, keeper, j, lim, x, y; 

loc p, q; 

char *ptr, *qptr; 

connections c, r; 

Boolean foo; 

for (len = 0; s->edge[len] != NULL; ++len) ; 
for (i = 0; i < len; ++i) { 
p = s->edge [i] ; 

foo = NEAR(p->x, 171,3) && NEAR(p->y, 182, 3); 
if (foo) PIF "%d,%d connected L=%x, R=%x\n", p->x, p->y, 

p->cleft & 255, p->cright & 255 ENDPIF 
for (goleft = 0, ptr = &p->cright; goleft <= 1; 

++goleft, ptr = &p->cleft) { 
keeper = f ind_range (*ptr, goleft); 
if (keeper < 0) 

continue; 
if (top == bottom) 

continue; 
++ callno; 
++ curcolor; 

pick_best_path (p, goleft, keeper, 0) ; 

} 

} 

} 

static void detect_edges (int **suspects, graphics_data *data) 
{ 

int i, j, *ptr, x, y, dir; 
loc p; 

for (j = 0; suspects [j] != NULL; ++ j ) { 
ptr = suspects [j]; 

for (i = 0; ptr[i] >= 0; i += 2) { 
X = ptr [i] ; 
y = ptr [i + 1] ; 

for (dir = 0; dir <= 1; ++dir) { 

p = f ollow_gradient (data, x, y, dir, 25); 
while (sp > 0) { 

p = stack [sp--] ; 

if (p->reached != 3) 

(void) f ollow_gradient (data, p->x, p->y/ 
2 - p->reached, 3) ; 

} • 

} 

} 



PIF "Xn" ENDPIFF 

} 

} 

Static loc pointbuffer [10000] ; 
static int log_2 [256] ; 

static consolidate (data) 
graphics_data *data; 

{ 

splotch s, n, lis; 

int i, inc, ix, con, dir, x, y, hi, lo; 
loc p, q; 

for (i = 1; i < 256; ++i) 

log_2 [i] = -2 ; 
log_2[0] = -1; 
for (i = 0; i < 8; ++i) 

log_2 [1 « i] = i; 
loop_over_list (s, data->f eatures) 

for (i = 0; s->edge[i] != NULL; ++i) 
s->edge [i] ->color = 0; 
curcolor = 1; 
lis = NULL; 

loop_over_list (s, data->f eatures) { 

for (i = 0; s->edge[i] != NULL; ++i) { 
p = s->edge [i] ; 
if (p->color >= curcolor) 
continue; 

if {p->cleft == 0 && p->cright == 0) 

continue; 
++ feature_id; 

PIF "Starting feature %d: %d,%d\n", feature_id, p->x, p->y 

ENDPIF 

for (inc = -1; inc <= 1; inc +=2) { 
ix = 5000; 

pointbuf f er [ix] = p; 
q = P; 

while (True) { 

q->color = curcolor + feature_id; 

q->f eature_id = feature_id; 

con = inc == -1 ? q->cleft : q->cright; 

dir = log_2 [con & 255] ; 

/* 

PIF " %d,%d: dir = (%x) %d\n", 

q->x, q->y, con & 255, dir ENDPIF 

*/ 

if (dir < 0) { 
if (dir < -1) 

PIF "Ambiguous point %d,%d: %c->%x\n", q->x, q- 

"RL" [inc == -1] , con & 255 ENDPIF 

break; 

} 

X - q->x + step_x[dir]; 
y = q->y - step_y[dir]; 
q = lookuploc(x, y) ; 



if (q->color >= curcolor) { 

PIF " Bailing at %d,%d due to color=%d\n", 
X, y, q->color ENDPIF 

break; 

} 

ix += inc; 

pointbuf fer [ix] = q; 

} 

*{iiic == -1 ? &lo : &hi) = ix; 

} 

n = new_splotch (&lis) ; 

n->edge = newarray (loc, hi - lo + 2) ; 

for (ix = lo; ix <= hi; ++ix) 

n->edge[ix - lo] = pointbuf fer [ix] ; 
n->edge[ix - lo] = NULL; 
q = pointbuf fer [lo] ; 
p = pointbuf fer [hi] ; 

PIF " Feature from %d(%d,%d) to %d (%d, %d) \n" , 
lo, q->x, q->y, hi, p->x, p->y ENDPIFF 

} 

} 

data->f eatures = lis; 

} 

static f ind_f olds (data) 
graphics_data *data; 

{ 

splotch s, r, orig; 

int i, len, dirl, dir2, dir3, angle, thl, th2, restart; 
loc p, q, t, u; 

loop_over_list (orig, data->f eatures) { 
s = orig; 

for (len = 0; s->edge[len] != NULL; ++len) ; 
if (len < edge_width * 3) 
continue; 

for (i = edge_width * 3; i < len; i += edge_width) { 
p = s->edge[i - edge_width * 3]; 
q = s->edge[i - edge_width * 2] ; 
t = s->edge[i - edge_width] ; 
u = s->edge [i] ; 

dirl = cartesian_to_j)olar (q->x - p->x, p->y - q->y) ; 
dir3 = cartesian_tojpolar (u->x - t->x, t->y - u->y) ; 
angle = (dir3 - dirl) & 7; 
if (absangle [angle] >= 3) { 

dir2 = cartesian_to_polar (t->x - q->x, q->y - t->y) ; 

thl = (dir2 - dirl) & 7; 

th2 = (dir3 - dir2) & 7; 

if (thl >= 5 && th2 >= 5) 

continue; // Backward fold --no sense breaking 

those 

PIF "Found a fold at %d,%d - %d,%d (dir %d -> %d -> %d)\n", 
q->x, q->y, 
t->x, t->y, 

dirl, dir2, dir3 ENDPIFF 



restart = i - edge_width; 



// Terminate old splotch early 

s->edge [restart - edge_width + 1] = NULL; 

// Should hunt for the best spot to break instead of just 
// snipping middle segment, but implement it later. 

// Start a new splotch after the break 
++ feature_id; 

r = new_splotch (&:data->f eatures) ; 
r->edge = &s->edge [restart] ; 

// Continue, looking for more folds in the new splotch 
s = r; 

len -= restart; 

i = edge_width * 2; 

} 

} 

} 

} 



typedef struct { 
splotch which; 
char dir; 
char len; 
loc p; 

} passthrurec, *passthru; 

typedef struct { 

passthrurec pt [2] ; 

short X, y; 
} phboxrec, *phbox; 

static put_injph (splotch s, phbox b, loc enters, loc leaves, int cnum) 
{ int dx, dy; int dir, i, angle, minangle, ix; passthru pt, ot; 

dx = leaves ->x - enters ->x; 

dy = leaves ->y - enters ->y; 

dir = cartesian_toj)olar (dx, - dy) ; 

i = dir / 2; 

if (cnum == ( (i - 1) & 3) ) 

return; 
if (cnum == ( (i - 2) & 3) ) 

return; 
minangle = 5; 
for (i = 0; i < 2; ++i) { 

pt = &b->pt [i] ; 

if (pt->which == s) 
angle = 0; 

else if (pt->len == 0) 
angle = 2; 

else 

angle = absangle [ (dir - pt->dir) & 7]; 
if (minangle > angle) { 
minangle = angle; 



ix = i; 

} 

} 

pt = &b->pt [ix] ; 
if (mag > pt->len) { 

PIF "Pigeoning %d,%d to %d,%d: s%d (%d ® %d [%d,%d])\n", 
enters->x, enters->y, leaves->x, leaves->y, 
s->id, mag, dir, dx, -dy ENDPIF 
pt->len = mag; 
pt->dir = dir; 
pt->which = s; 
pt->p = enters; 

} 

} 

static splotch chase_mesas (splotch s) 

{ 

splotch p; 

for (p = s; p ! 
if (p->mesa = 
break; 
return (p) ; 

} 

static group (si, s2) 
splotch si, s2; 

{ 

splotch r, p; 

r = chase_mesas (si) ; 
p = chase_mesas (s2) ; 

PIF "Grouping %d(%d) with %d(%d)\n", sl->id, r->id, s2->id, p->id 
ENDPIF 

if (p->id < r->id) 

r->mesa = p; 
else 

p->mesa = r; 

} 

static record_pairing (phbox box) 

{ 

pairing p; 
splotch si, s2; 

si = box->pt [0] .which; 
s2 = box- >pt [1] .which; 
if (sl->id > s2->id) { 

si = box- >pt [1] .which; 

s2 = box->pt [0] .which; 

} 

loop_over_list (p, sl->partners) { 
if (p->sl == si ScSc p->s2 == s2) 
break; 

} 

if (p == NULL) { 

p = new{pairingrec) ; 



= NULL; p = p->mesa) 
= NULL I I p->mesa == p) 



p->next = si- >partners; 
sl->partners = p; 
p->Sl = si; 
p->s2 = s2; 
p->COunt = 0; 

} 

++ p->count; 
if (p->count == 2) 
group (si, s2) ; 

} 

static pairing last_partner (pairing p) 

{ 

for ( ; p->next != NULL; p = p->next) ; 
return (p) ; 

} 

static bipartition (splotch s) 
{ 

pairing p, q; 

int max, pi, p2, newpl, newp2; 
Boolean useful; 



q = NULL; 
max = 0; 

loop_over_list (p, s->partners) 
if (maix < p->count) { 
max = p->count; 
q = P; 

} 

if (max < 2) 

return; 
q->sl->polarity = max; 
q->s2->polarity = - max; 
useful = True; 
while (useful) { 

useful = False; 

loop_over_list (q, s->partners) { 
newpl = newp2 = 0; 
pi = q->sl->polarity; 
p2 = q->s2->polarity; 
if (pi == 0) { 

newpl s= - p2; 
} else if (p2 == 0) { 

newp2 = - pi; 
} else if ((pi > 0) == (p2 > 0) ) { 

if {(pi > 0) ? pi > p2 : pi < p2) 
newp2 - - pi; 

' else 

newpl = - p2 ; 
} else { 

if {(pi > 0) ? pi > - p2 : pi < - p2) 
newp2 = - pi; 

else 

newpl = - p2; 

} 

if (newpl > q->count) newpl = q->count; 



if (newpl < - q->count) newpl = - q->count; 
if (newp2 > q->count) newp2 = q->count; 
if (newp2 < - q->count) newp2 = - q->count; 
if (abs(pl) < abs (newpl)) { 

q->sl->polarity = newpl; 

useful = True; 

} 

if (abs{p2) < abs{newp2)) { 
q->s2->polarity = newp2; 
useful = True; 

} 

} 

} 

} 

static pigeon_hole (data) 
graphics_data *data; 

{ 

splotch s; 

int nx, ny, i, x, y, ix, iy, j, k, id, cnum; 
phbox boxes, box, corner [5]; 
int stepx[5], stepy[5] ; 
loc p, entry [5] ; 
pass thru pt; 

PIF "Pigeonholing\n" ENDPIFF 
nx = (hi_x " lojx. - 1) / box_size + 1; 
ny = (hi_y - lo_y - 1) / box_size + 1; 
boxes = newarray (phboxrec, nx * ny) ; 
for (i = nx * ny; --i >= 0; ) { 
box = &boxes[i]; 

box->pt [0] -len = box->pt [1] . len = 0; 
box->x = i % nx * box_size + lo_x; 
box->y = i / nx * box_size + lo_y; 

} 

for (i = 0; i < 4; ++i) { 

stepx[i] = step_x[i * 2 + 1] * edge_width; 
stepy[i] = step_y[i * 2 + 1] * edge_width; 

} 

stepx[4] = stepy[4] = 0; 
loop_over_list (s, data->f eatures) { 
id = s->id; 

for (cnum = 0; cnum < 5; corner [cnum++] = NULL); 
for (i = 0; s->edge[i] != NULL; ++i) { 
p = s->edge [i] ; 

for (cnum = 0; cnum < 5; cnum++) { 

X = (p->x - lo_x + stepx[cnum]) / box_size; 
y = (P'>y - lo_y - stepy[cnum]) / box_size; 
if (x >= nx I I y >= ny) 
box = NULL; 

else 

box = &boxes [y * nx + x] ; 
if (box != corner [cnum] ) { 

if (corner [cnum] != NULL) { 

put_in_j)h (s, corner [cnum] , entry[cnum], p, cnum); 

} 

corner [cnum] = box; 



} 



entry [cnum] = p; 

} 

} 

} 

for (cnum = 0; cnum < 5; cnum++) { 
if (corner [cnum] 1= NULL) { 

put_in__ph(s, corner tcnum] , entry [cnum], p, cnum); 

} ^ 

} 

for (y = 0/ y < ny; ++y) { 
for (x = 0; X < nx; ++x) { 

box = &boxes [y * nx + x] ; 
if (box->pt [0] .len == 0) 

continue; 
PIF "%d,%d - %d,%d:", 

x * box_si2e + lo_x, y * box_size + lo_y, 
(x + 1) * box_size + lo_x, (y + 1) * box_size + lo_ 
ENDPIF 

for (i = 0; i < 2; ++i) { 
pt = &box->pt [i] ; 
if (pt->len == 0) 
break; 

PIF " (s%d dir=%d, len=%d)", 

pt->wliicli->id, pt->dir, pt->len ENDPIF 

PIF "\n" ENDPIF 
if (i == 2) 

record_j)airing (box) ; 

} ^ 

loop_over_list (s, data->f eatures) 

s->mesa = cliase_mesas (s) ; 
loop_over__list (s, data- >f eatures) 
if (s->mesa != s) 

if (s->partners != NULL) { 

last_partner (s->partners) ->next = s->mesa->partners 
s->mesa->partners = s->partners; 
s->partners = NULL; 

} 

loop_over_list (s, data->f eatures) 
if (s->mesa == s) 
bipartition(s) ; 
loop_over__list (s, data->f eatures) { 

PIF "Splotch %d is in mesa %d (pol=%d)\n", 

s->id, s->mesa->id, s->polarity ENDPIF 



static orient_mesa (splotcli *pile) 

{ 

int i, j, k, pos, plus_inside, minus_inside, cutoff; 
int *min, *max; *w]iich; 
splotch s; 

int plus_bbox [4] , minus_bbox [4] , dif f [4] ; 
loc p; 



for (i = 0, which = plus_bbox; i < 2; which = minus_bbox) { 

min = which; 
max = which + 2 ; 
min[0] = min[l] = 1000000; 
max[0] = max[l] = -1; 

} 

for (i = 0; pileti] != NULL; ++i) { 
s = pile [i] ; 

which = s->polarity > 0 ? plus_bbox : minus_bbox; 
min = which; 
max = which + 2 ; 

for (j = 0; s->edge[j] != NULL; ++j) { 
p = s->edge [j] ; 

for (k = 0, pos = p->x; k < 2; ++k, pos = p->y) { 
if {min[k] > pos) 

mintk] = pos; 
if (max[k] < pos) 

max[k] = pos; 

} 

} 

} 

for (i = 0; i < 4; ++i) 

dif f [i] = plus_bbox[i] - minus_bbox [i] ; 
for (i = 2; i < 4; ++i) 

diff [i] = - diff [i] ; 

cutoff = 2; 

plus_inside = minus_inside = 0; 
for (i = 0; i < 4; ++i) 
if {diff[i] < - cutoff) 

++ minus_inside; 
else if (diff [i] > cutoff) 
++ plus_inside; 

which = plus_inside < minus_inside ? plus__bbox : minus_bbox; 
s = pile [0] ->mesa; 

PIF "Bounding box of mesa %d is %d,%d to %d,%d\n", s->id, 

which[0], which[l], which[2], which[3] ENDPIF 
PIF " Red inside: %d, Blue inside: %d\n", 

plus_inside, minus_inside ENDPIFF 
if {minus_inside > plus_inside) { 
PIF " Reversing\n" ENDPIFF 
for (i = 0; pile [i] 1= NULL; ++i) { 
s = pile [i] ; 

s->polarity = - s->polarity; 

} 

} 

} 

static decide_inside_outside (data) 
graphics_data *data; 

{ 

splotch s, t; 

int count; 
splotch pile [100] ; 

loop_over_list (s, data->f eatures) 
if {s->mesa == s) { 
count = 0; 



loop_over_list (t, data->features) 
if (t->mesa == s) 

if (t->polarity <= -2 | | t->polarity >= 2) 
pile [count++] = t; 
pile [count] = NULL; 
if (count >= 2) 

orient_mesa (pile) ; 

} 

} 

exatnine_cross sections (graphics data *data) 

{ 

windo sem; 

int xsize, ysize, count, i, j; 
double f, separation; 
splotch s; 
int *suspects [50] ; 



data->f eatures = NULL; 
if ( (int) cfp == -1) 

cfp = fopen("cross_section" , "w"); 

theDisp=XtDisplay (data->seTn->thiscanvas) ; 

data- >f eatures = NULL; 

init_hash_table 0 ; 

identify_border_areas (data) ; 

sem = data->sem; 

xsize = hi__x - lo_x; 

ysize = hi_y - lo_y; 

separation = 1.0 / 6; 

for (i = 0; i <= HIST_SIZE; histogram [i++] = 0); 
count = 0; 

for (f = separation * 0.5; f < 0.99; f += separation) 

suspects [count++] = cross_section(data, lo_y + (int) (ysize * f ) , 

1); 

for (f = separation * 0.5; f < 0.99; f += separation) 

suspects [count++] = cross_sect ion (data, lo_x + (int) (xsize * f), 

0); 

suspects [count] = NULL; 
analyze_pair_statistics () ; 
detect_edges (suspects, data) ; 
loop_over__list (s, data- >f eatures) 

break_links (s) ; 
consolidate (data) ; 
f ind_f olds (data) ; 
pigeon_hole (data) ; 
decide_inside_outside (data) ; 
PIF "\n" ENDPIFF 
/* 

printf ("hcalls = %d, hpoints = %d, hsteps = %d, hmax = %d\n", 

hcalls, hpoints, hsteps, hmax); 
printf ("hrepeats = %d\n", hrepeats) ; 
*/ 

} 



Subclaim 1-e excerpt from callbacks. c May 25, 1999 - 



This file contains a procedure for drawing the contours of an 
alim image (an aerial image or latent image) onto a digitized 
wafer pattern image. When both the alim contours and the 
detected pattern edges are drawn onto the same pattern image, 
they are overlaid on each other. 

/* 

% Various Callback Routines and Drawing Functions for SmartSEM 
% 

% Dec. 1995 

*/ 



#include 
#include 
ftinclude 
ttinclude 



<stdio.h> 

<math.h> 
"Screens .h" 
"plot .h" 



extern int DrawFillMaskFlag, DrawGrid; 

extern MaskRegion mrd; 

extern int npregions; 

extern int DrawImagelnMaskFlag; 

extern Pixel red_j>ixel; 

extern int 

color_map_type, color_mask, color_label, color_SEM_contour ; 
extern int width_label , width_SEM_contour , width_mask; 

extern double contour_value; 
extern double contour_list [] ; 

extern windo newwindo ( ) ; 



Widget trackerP, 
trackerW, 
valueDisplay, 
trackerD; 



static 
static 

int ixoffset, iyoffset, y_top; 
double x_f actor, y_factor/ 



/* widgets for tracking cursor 
/* pixel, world x,y and */ 
/* function min and max values 

/* function value at point*/ 



graphics_data *contour_data; 
Widget contour__widget; 



^•k1t*it1t****ii*-k-k-kitir1c-k****ic 

#include <sys/time.h> 
hrtime_t start, end; 



Contour Drawing Functions 



static int pixcount; 

static int matrix_size = 0; 

static float *2_matrix; 

static Display *contour_dpy ; 

static float xavg, yavg, zavg; 



static Boolean doSIM; 
static float *x, *y; 



#define FAREF(i,j) ((float) z__matrix [ ( j ) *nx+ (i) ] ) 

world_to_pixel_conversion (graphics data *data, Boolean sem) 
{ 

windo w; 

if (sem) { 

w = data->sem; 

y_top = w->yres; 
} else { 

w = data->aim; 

y_top = data->height; 

} 

x_f actor = 1.0 / w->pixel.x; 
y_f actor = 1.0 / w->pixel.y; 

ixoffset = (int) (w->ll.x * x_factor + 0.5); 
iyoffset = - (int) (w->ll.y * y factor + 0.5); 

} 

/* Draw a line between two points, given pixel coordinates. */ 
static void contour_seg (int ix2, int iy2, int ix, int iy) { 
XDrawLine (contour_dpy , XtWindow (contour_widget ) , 
cont our_dat a - >gcSEMContour , 
ix2,iy2,ix,iy) ; 

} 

/* Draw a line between two points, given world coordinates. */ 
static connect joints (xl, yl, x2, y2) double xl, yl, x2, y2; { int ix, 
iy, ix2, iy2; 



ix = (int) (xl * x_f actor +0.5) - ixoffset; 

iy = y_top - (int) (yl * y_f actor +0.5) - iyoffset; 

ix2 = (int) (x2 * x_f actor +0.5) - ixoffset; 

iy2 = y_top - (int) (y2 * y_factor +0.5) - iyoffset; 

contour_seg(ix2, iy2, ix, iy) ; 

} 

/* Draw the line at a given z-value across a triangle. */ static 
int plotContourSegment (float level, float xl, float yl, float zl, 
^ float x2, float y2, float z2) 

float x3, y3, z3; 
graphics_data *data ; 
float X, y, factor, temp; 
int ix , iy , ix2 , iy2 ; 
float xt,yt; 

data = contour_data; 
x3 = xavg; 
y3 = yavg; 
z3 = zavg; 



if (21 >= 22) { 

temp = xl; xl = x2; x2 = temp; 
temp = yl; yl = y2; y2 = temp; 
temp = zl; zl = z2; z2 = temp; 

} 

if(z2 > z3) { 
if(z3 < zl) { 

temp = xl; xl = x3 ; x3 = temp; 
temp = yl; yl = y3; y3 = temp; 
temp = zl; zl = z3; z3 = temp; 

temp = x2; x2 = x3; x3 = temp; 
temp = y2; y2 = y3 ; y3 = temp; 
temp = z2; z2 = z3; z3 = temp; 

} 

/* zl <= 22 <= z3 */ 
if (level < zl I I level > z3) { 
return -1; 

} 

if ((level == zl) && (zl == z2) && (z2 != z3) ) { 
/* draw a line from xl, yl to x2, y2 */ 
connect^points (xl , yl , x2 , y2 ) ; 
return 1; 

} 

if ((level == z3) && (z3 == z2) && (zl != z2)) { 
/* draw a line from x2, y2 to x3, y3 */ 
connect_points (x3 , y3 , x2 , y2 ) ; 
return 1; 

} 

if (level ==22) { 

/* we know that z2 is definitely between zl and z3 

draw a line from P2 to point on segment between PI and P3 

*/ 

if (23 == zl) 

return 0; 
factor = (level - zl)/(z3-2l); 
X = xl + (x3-xl) * factor; 
y = yl + (y3-yl) * factor; 
connect_points (x, y, x2, y2) ; 
return 1; 

} 

if (level < z2) { 

if (z2 == zl II z3 == zl) 

return 0; 
/* P2 and P3 are above, PI below */ 
factor = (level - zl) / (z2 - zl) ; 
xt = xl + (x2 - xl) * factor; 
yt = yl + (y2 - yl) * factor; 

factor = (level - zl) / (z3 - zl) ; 
X = xl + (x3 - xl) * factor; 
y = yl + (y3 - yl) * factor; 
connectjpoints (x, y, xt, yt) ; 
} else { 



if (z3 == zl II z3 == z2) 
return 0; 

/* PI and P2 are below, P3 above */ 
factor = (level - zl) / (z3 - zl) ; 
xt = xl + (x3 - xl) * factor; 
yt = yl + (y3 - yl) * factor; 

factor = (level - z2) / {z3 - z2) ; 
X = x2 + (x3 - x2) * factor; 
y = y2 + {y3 - y2) * factor; 
connect_points (x, y, xt, yt) ; 

return 1; 

} 

/* This system draws a contour line over a display square by 
identifying which sides of the square cross the contour 
line. If zero do, there's no line to draw. If all four 
do, this method doesn't work, so Bob Pack's old pyramid 
method is used. But if two sides cross the line, this 
technique identifies them by table- lookup based on which 
corners are above and below the contour level, and then 
draws a line between the two crossing points, pgf, 4/99 */ 

#define NO__CONTOUR 5, 5, 5, 5 
#define PYRAMID 4, 4, 4, 4 

typedef struct { 

int lol, hil, lo2, hi2; 
} edgepairrec; 

static edgepairrec edgepairs [17] = { // Corner Numbers 

// Corner: 3210 



N0_CONTOUR, 


II 


LLLL 


1, 0, 2, 


0, 


II 


LLLH 


0, 1, 3, 


1, 


II 


LLHL 


2, 0, 3, 


1, 


II 


LLHH 


0, 2, 3. 


2, 


II 


LHLL 


3, 2, 1, 


0, 


II 


LHLH 


PYRAMID, 




II 


LHHL 


3, 2, 3, 


1, 


II 


LHHH 


1, 3, 2, 


3, 


II 


HLLL 


PYRAMID, 




II 


HLLH 


0, 1, 2, 


3, 


II 


HLHL 


2, 0, 2, 


3, 


II 


HLHH 


0, 2, 1, 


3, 


II 


HHLL 


1, 0, 1, 


3, 


II 


HHLH 


0, 1, 0, 


2, 


II 


HHHL 


NO_CONTOUR, 


II 


HHHH 



-1}; 

/* Find the points where a contour level crosses two edges of 
a display square, and draw a line between them. */ 
static connect_edges (edgepairrec *ep, int i, int j, double level) { int 
ixlo, ixhi, iylo, iyhi; double xlo, xhi, ylo, yhi, zlo, zhi, frac; 
double xl, yl, x2, y2; 



ixlo = i + (ep->lol Sc 1) ; 
iylo = j + {ep->lol >> 1) ; 
ixhi = i + {ep->hil & 1) ; 
iyhi = j + {ep->hil >> 1) ; 
xlo = X [ixlo] ; 
xhi = x tixhi] ; 
ylo = y [iylo] ; 
yhi = y [iyhi] ; 
zlo = FAREF(ixlo, iylo) / 
zhi = FAREF(ixhi, iyhi) ; 
frac = (level - zlo) / (zhi - zlo) ; 
if (frac > 1.0 II frac < 0.0) 
return; 

xl = xlo + frac * (xhi - xlo) ; 
yl = ylo + frac * (yhi - ylo) ; 

ixlo = i + {ep->lo2 & 1) ; 
iylo = j + (ep->lo2 >> 1) ; 
ixhi = i + (ep->hi2 & 1) ; 
iyhi = j + {ep->hi2 >> 1) ; 
xlo = X [ixlo] ; 
xhi = X [ixhi] ; 
ylo = y [iylo] / 
yhi = y [iyhi] ; 
zlo = FAREF(ixlo, iylo) ; 
zhi = FAREF ( ixhi , iyhi ) ; 
frac = (level - zlo) / (zhi - zlo) ; 
if (frac > 1.0 || frac < 0.0) 
return; 

x2 = xlo + frac * (xhi - xlo) ; 
y2 = ylo + frac * (yhi - ylo) ; 

connectjpoints (xl , y 1 , x2 , y2 ) ; 

} 

/* Draw the segments for a set of contour lines that pass through 
a one-grid-wide vertical stripe of an AIM or SEM. */ 
static draw_contour_set (double *contourset, int i. Boolean drawgrid) { 
Pointjptr data_j>tr; int j, k; 
float level, zij , zipj , zijp, zipjp; 

double contour2 [2] ; 

int top_above [100] , bottom_above [100] ; 
int shape; 
edgepairrec *ep; 
windo w; 

for (j = 0; j < ny - 1; j++) 

if (y[j+l] > contour_data->aim->ll.y) 
break; 
zijp = FAREF (i, j ) ; 
zipjp = FAREF{i+l, j); 

for ()c = 0; contourset [}c] > 0.0; ++)c) { 
level = contourset [k] ; 

top_above[k] = (zijp > level) | ((zipjp > level) << 1); 

} 

for ( ; j < ny - 1; j++) { 

if (y[j] >= contour_data->aiTn->ur .y) 



break; 
zij = zijp; 
zipj = zipjp; 
zijp = FAREF{i, j+1); 
zipjp = FAREF(i+l, j+1) ; 



yavg = (y [j] + y [j+1] ) /2 .0; 
zavg = (zij+zipj+zijp+zipjp) /4 . 0; 
for (k = 0; contourset [k] > 0.0; ++k) { 
level = contourset [k] ; 

/* Determine which corners of the current display square 
are above the contour level and which are below. */ 
bottom_above [k] = top_above [k] ; 

top_above [k] = (zijp > level) | ((zipjp > level) << 1); 
shape = bottom_above [k] | (top_above [k] << 2); 
ep = &edgepairs [shape] ; 



if (ep->lol == 5) 
continue; 

if (ep->lol == 4 I I pixcount > 15) { 

/* Pyramid method. Chop the display square into four 

triangles and compute the contour line for each. */ 
plotContourSegment { level , 

x[i], y[j], zij, x[i], y [j+1] , zijp) ; 
plotContourSegment (level, 

x[i], y[j], zij, x[i+l],y[j], zipj); 
plotContourSegment (level, 

x[i+l],y[j], zipj, x[i+l] ,y [j+1] , zipjp) ; 
plotContourSegment (level, 

x[i], y[j+l],zijp, x[i+l] ,y[j+l] , zipjp) ; 

} else { 

connect_edges (ep, i, j, level); 

} 

} 

if (drawgrid) { 
if (doSIM) 

w = contour_data->aim; 
else 

w = contour_data->sem; 
connect_points(x[i+l] - w->pixel.y, y[j+l], x[i+l], y[j+l]); 

} 



/* 

* Draw Contour Lines RC Pack Nov 1997 
*/ 

void drawContours (Widget w, graphics^data *data, int callFunction) { 

Point jptr data_j)tr; 

int i,ix,iy; 

float xlas t , ylast ; 

Boolean generatel, generate2; 

double contour2 [2] ; 



if {data->image == NULL | | w == NULL) 

return; 
doSIM = callFunction != 1; 



if (doSIM) { 

// generatel = data->drawContourOnSIMl/ 
// generate2 = data->drawContourOnSIM2 ; 

generatel = XmToggleButtonGetState (data->aim->dependents [1] ) ; 
generate2 = XmToggleButtonGetState (data->aiTn->dependents [3] ) ; 

} else { 

generatel = XmToggleButtonGetState (data->sem~>dependents [1] ) ; 
generate2 = XmToggleButtonGetState (data->sem->dependents [2] ) ; 

if (! generatel && ! generate2) 
return; 

switch_image (data, 5) ; 

contour_data = data; 
contour_widget = w; 
contour_dpy = XtDisplay{w) ; 
if (generatel) 

parse_contour_list (data, generate2) ; 

if (XtlsRealized(w) ) { 
// llApr98 

if { width_SEM_contour== - 1 ) 

width_SEM_contour=0; // A kludge until initalized 
XSetLineAttributes (contour_dpy, data->gcSEMContour, 

(unsigned) width_SEM_contour , LineSolid, CapButt , JoinMiter ) ; 



X = malloc (sizeof (float) * nx) ; 
y = malloc (sizeof (float ) * ny) ; 
if (nx * ny > matrix_size) { 
matrix_size = nx * ny; 

z_matrix = malloc (sizeof (float) * matrix size); 

} 

xlast = ylast = -9999.0; 
ix=iy=0; 

datajptr = data->image->aimvalues; 
for (i=0; i<nx*ny ; i++) { 
z_matrix [i] =data_ptr->z; 
if (data_ptr->x > xlast) 

xlast=x [ix++] =datajptr->x; 
if (data_ptr->y > ylast) 

ylast=y [iy++] =dataj)tr->y; 
data_ptr++; 

} 

contour2[0] = data->contour_value; 
contour2 [1] = -1.0; 

if (doSIM) 

pixcount = data->image->del .y / contour_data->aim->pixel .y; 
else 

pixcount = data -> image ->del.y / contour_data->sem->pixel .y; 
world_tojpixel_conversion{data, ! doSIM) ; 



// Contour Plot 
for (i = 0; i < nx - 1; i++) 
if {x[i+l] > data->aim->ll .x) 
break; 
start = gethrtime {) ; 
for ( ; i < nx - 1; i++) { 
if (x[i] >= data->aim->ur .x) 
break; 

xavg = (x[i] + x[i+l]) / 2.0; 

if (generatel) { 

XSet Foreground (contour_dpy , data->gcSEMContour, 
color_SEM__contour) ; 

draw_contour set (contour list, i, DrawGrid) ; 

} 

if (generate2) { 

XSe t Foreground (contour_dpy, data->gcSEMContour , red_pixel) ; 
draw_contour_set (contour2, i, DrawGrid && ! generatel); 

} 

end = gethrtime {) ; 
// printf ("Display time = %.2f msec\n", (end - start) * l.Oe-6); 
} 

Subclaim 1-f excerpt from semedge.c June 1, 1999 



The following addition to semedge.c measures and reports the 
average 

and root -mean- square distances from the contours in an alim image 
to the detected edges . 

/************************ SEM/AIM Comparison Functions 

ititirliicffkftitleitiriclrlt********^ 

Static transf orm_to_orig (xl, yl, px, py, data) 
double xl, yl; 
int *px, *py; 
graphics_data *data; 

{ 

double X, y; 
int lev; 



X = (xl * x_factor) - (ixof fset + sem__of f set) ; 

y = y_top - (iyoff set + sem_of fset) - (yl * y_f actor) ; 

if (data->f lipX) 

X = data->sem->xres - x; 
if (data->f lipY) 

y = data->sem->yres - y; 
lev = data->curpm->zoom_level; 
x = X * ZOOM_STEPS / lev; 
y = y * ZOOM_STEPS / lev; 
*px = (int) (x + 0.5) ; 
*py = (int) (y + 0.5) ; 



static find_nearby_snakes (float *Beg, graphics_data *data) 

double sxl, syl, sx2, sy2 ; 
int xl, yl, x2, y2 ; 

int ix, iy, mx, my, x, y, i, j, count, length; 

int lox, loy, hix, hiy; 

int nx, ny; 

int minmagl, minmag2; 

phbox box; 

snake s; 

loc p; 

int distance, squaresum, metric; 



sxl = 


seg [0] ; 


syl = 


seg[l] ; 


sx2 = 


seg [2] ; 


sy2 = 


seg [3] ; 



transform_to_orig(sxl, syl, ScXl, &yl, data); 
transf orm_to_orig ( sx2 , sy2 , &x2 , &y2 , data) ; 

if (xl == x2 ScSc yl == y2) 
return; 

cartesian_to_polar (xl - x2, yl - y2) ; 
length = mag; 

nx = (hi_x - lo_x - 1) / box_size + 1; 

ny = {hi_y - lo_y - 1) / box_size + 1; 

mx = (xl + x2) >> 1; 

my = (yl + y2) >> 1; 

ix = (mx - lo_x) / box_size; 

iy = (my - lo_y) / box_si2e; 

lox = (ix - RANGE) * box_size + lo_x; 

loy = (iy - RANGE) * box_size + lo_y; 

hix = lox + (2 * RANGE + 1) * box_size; 

hiy = loy + (2 * RANGE + 1) * box_size; 

PIF "Segment %g,%g - %g,%g uM (Sem-edge database coords: %d,%d - 
%d,%d)\n", 

sxl, syl, sx2, sy2, 
xl, yl, x2, y2 
ENDPIF 

count = 0; 

for (x = ix - RANGE; x <= ix + RANGE; ++x) { 
if (x < 0 I I X >= nx) 
continue; 

for (y = iy - RANGE; y <= iy + RANGE; ++y) { 
if (y < 0 I I y >= ny) 

continue; 
box = &:boxes [y * nx + x] ; 
for (i = 0; i < 2; +-i-i) { 
if (box->pt [i] .len == 0) 

break; 
s = box- >pt [i] .which; 
if (s->dir != 0 I I s->polarity > -2) 

continue; 
semedge [s->dir = ++count] = s; 

} 



} 

} 

lox = (xl < x2 ? xl : x2) - box_size; 

hix = (xl > x2 ? xl : x2) + box_size; 

loy = (yl < y2 ? yl : y2) - box_size; 

hiy = (yl > y2 ? yl : y2) + box_size; 

distance = 1000; 

squaresum = 1000000; 

for (i = 1; i <= count; ++i) { 

s = semedge [i] ; 

s->dir = 0; 

minmagl = minmag2 = 1000000; 
for (j = 0; s->edge[j] != NULL; ++ j ) { 
p = s->edge [j] ; 

if (p->x < lox II p->x > hix II p->y < loy | | p->y > hiy) 
continue; 

cartesian_tojpolar (p->x - xl, p->y - yl) ; 
if (mag < minmagl) 

minmagl = mag; 
cartesian__tojpolar (p->x - x2, p->y - y2) ; 
if (mag < minmag2) 

minmag2 = mag; 

} 

PIF " distance to s%d: %d %d\n", s->id, minmagl, minmag2 

ENDPIF 

metric = minmag2 + minmagl; 
if (distance > metric) 
distance = metric; 
metric = {minmag2 - minmagl) ; 

metric = metric * metric + minmag2 * minmagl * 3; 
if (squaresum > metric) 
squaresum = metric; 

} 

if (distance > 100) 

return; 
totallength += length; 
totalarea += length * distance; 
totalsquares += length * squaresum; 



static void measure_edge_contour_dif f erence (graphics_data *data) { 
snake s; int i, segments; float *seg; double gap, pixsize; 

if (data->features->how_close >= 0.0) // if already computed 
return; 

boxes = data->f eatures->boxes; 

printf ( "\nComputing the gap between the edges and the contour 
lines . . . \n" ) ; 

for (i = 0; contour_store [i] > -l.e6 || 

contour_store [i + 1] > -1.0e6; i += 2) { 

segments = i / 4; 

totallength = totalarea = totalsquares = 0; 
for (i = 0; i < segments; ++i) { 

seg = &contour_store [i * 4] ; 

f ind_nearby_snakes (seg, data) ; 

} 



pixsize = {data->sem->pixel.x + data->sem->pixel .y) / 2; 
totalarea >>= 1; 
totalsquares /= 3/ 

gap = sqrt { (double) totalsquares / totallength) ; 
gap = gap * data->curpm->zoom_level / ZOOM_STEPS; 
data->f eatures->how_close = gap * pixsize * 1000; 
printf(" ...RMS separation = %.lf nM (%.2f pixels)\n", 
data->f eatures->how_close, gap) ; 

gap = (double) totalarea / totallength; 

gap = gap * data->curpm->zoom_level / ZOOM_STEPS; 

PIP "Average separation = %-lf nM (%.2f pixels) \n", 

gap * pixsize * 1000, gap ENDPIF 
printfC ...Average separation = %.lf nM (%.2f pixels)\n", 
gap * pixsize * 1000, gap); 

} 

report_edge_contour_difference (graphics data *data) 

{ 

char report [80] ; 

measure_edge_contour_dif f erence (data) ; 

sprintf (report, "RMS separation = %.lf nM", data->f eatures- 
>how_close) ; 

show_i f ^changed ( " C " , report ) ; 

} 

clear_how_close {graphics_data *data) 

{ 

if (data->f eatures == NULL) 
return; 

data- >f eatures ->how_close = -1.0; 

} 

Subclaim 1-g -- excerpt from semedge.c Feb. 29, 2000 



In this and subsequent versions, measure_edge_contour_dif f erence ( ) 
prints Molotof 's best guess at what the next iteration's value of 
the alim image contour should be. The contour is the scaled 
reciprocal of a mask processing parameter the amount of light 
to expose the wafer with. 

static int f_cartesian_tojpolar (dx, dy) 

double dx, dy; 

{ 

static int octant_map [8] = {O, 1, 3, 2, 7, 6, 4, 5}; 
double adx, ady; 
int ix; 

if (dx < 0) 

ix = 2, adx = - dx; 
else 

ix = 0, adx = dx; 
if (dy < 0) 

ix += 4, ady = - dy; 
else 



ady = dy; 
if (adx < ady) 

++ ix, adx *= 0.5; 
else 

ady *= 0.5; 
fmag = adx + ady; 
return {octant_map [ix] ) ; 

} 

static transf orm_f rom_orig (xl, yl, px, py, data) 
double xl, yl; 
double *px, *py; 
graphics_data *data; 

{ 

double X, y; 
int lev; 

lev = data->curpm->zoom_level; 
X = xl * lev / (double) ZOOM_STEPS; 
y = yl * lev / (double) ZOOM_STEPS; 
if (data->flipX) 

X = data->sem->xres - x; 
if (data->flipY) 

y = data->sem->yres - y; 
X += (ixoffset + sem_of fset) ; 
y -= y_top - (iyof f set + sem_of f set) ; 
X /= x_factor; 
y /= - y_factor; 
*px = x; 

*py = y; 

} 

static find_nearby_snakes (float *seg, graphics_data *data) 

{ 

double sxl, syl, sx2, sy2; 

double xl, yl, x2, y2 , mx, my, length; 

int ix, iy, x, y, i, j, count, dir; 

int lox, loy, hix, hiy; 

int nx, ny; 

double minmagl, minmag2; 
phbox box; 
snake s; 

loc p, pi, p2, cpl, cp2; 

double distance, squaresum, metric, z, slopel, slope2, pull_in, 
push_out; double dzdx, dzdy, grad; 

8x1 = seg [0] ; 
syl = seg [1] ; 
sx2 = seg [2] ; 
sy2 = seg [3] ; 

transform_to_orig (sxl, syl, &xl, &yl, data); 
transf orm_to_orig (sx2 , sy2, &x2, &y2, data); 
dir = f_cartesian_tojpolar (xl - x2, yl - y2); 
length = fmag; 

if (length < 0.001) // length is in pixels, 
return; 



nx = (hi_x - lo_x - 1) / box__size + 1; 
ny = (hi_y - lo_y - 1) / box^size + 1; 
mx = (xl + x2) * 0.5; 
my = (yl + y2) * 0.5; 
ix = (int) ( (mx - lo_x) / box__size) ; 
iy = (int) ({my - lo_y) / box_size) ; 
lox = (ix - RANGE) * box_size + lo_x; 
loy = (iy - RANGE) * box_size + lo_y; 
hix = lox + (2 * RANGE + 1) * box_size; 
hiy = loy + (2 * RANGE + 1) * box_size; 

PIF "Segment %g,%g - %g,%g uM (Sem-edge database coords: %g,%g - 
%g,%g)\n", 

sxl, syl, sx2, sy2, 
xl, yl, x2, y2 
ENDPIF 

count = 0; 

for (x = ix - RANGE; x <= ix + RANGE; ++x) { 
if (x < 0 I I X >= nx) 
continue; 

for (y = iy - RANGE; y <= iy + RANGE; ++y) { 
if (y < 0 I I y >= ny) 

continue; 
box = &boxes [y * nx + x] ; 
for (i: = 0; i < 2; ++i) { 
if (box->pt [i] . len == 0) 

break; 
s = box->pt [i] .which; 
if (s->dir != 0 | | s->polarity > -2) 

continue; 
semedge [s->dir = ++count] = s; 

} 

} 

} 

lox = (int) (xl < x2 ? xl : x2) - box_size; 
hix = (int) (xl > x2 ? xl : x2) + box_size; 
loy = (int) (yl < y2 ? yl : y2) - box_size; 
hiy = (int) (yl > y2 ? yl : y2) + box_size; 
distance = 1000; 
squaresum = 1000000; 
pi = p2 = cpl = cp2 = NULL; 
for (i = 1; i <= count; ++i) { 

s = semedge [i] ; 

s->dir = 0; 

minmagl = minmag2 = 1000000.0; 
for (j = 0; s->edge[j] 1= NULL; ++j ) { 
p = s->edge [j] ; 

if (p->x < lox I I p->x > hix I I p->y < loy | | p->y > hiy) 
continue; 

f__cartesian_to_polar (p->x - xl, p->y - yl) ; 
if (fmag < minmagl) { 

minmagl = fmag; 

cpl = p; 

} 

f_cartesian_to_polar (p->x - x2, p->y - y2) ; 
if (fmag < minmag2) { 

minmag2 = fmag; 

cp2 = p; 



} 

} 

PIF " distance to s%d: %d %d\n", s->id, minmagl, minmag2 

ENDPIF 

metric = minmag2 + minmagl; 
if (metric < distance) 
distance = metric- 
metric = {minmag2 - minmagl) ; 

metric = metric * metric + minmag2 * minmagl * 3; 

if (metric < squaresum) { 
squaresum = metric- 
pi = cpl; 
p2 = cp2; 

} 

} 

if (distance > 100) 

return; 

accumulate_vector (xl, yl, pl->x - xl, pl->y - yl, length, dir) 
accumulate_vector (x2 , y2, p2->x - x2, p2->y - y2, length, dir) 
total length += length; 
totalarea += length * distance; 
totalsquares += length * squaresum; 

f_cartesian_to_polar (pl->x - xl, pl->y - yl) ; 
minmagl - fmag; 

f_cartesian_to_polar (p2->x - x2, p2->y - y2) ; 
minmag2 = fmag; 

transform_f rom_orig ( (double) pl->x, (double) pl->y, &xl, Scyl, 
data) ; 

transform_from_orig{ (double) p2->x, (double) p2->y, &x2, &y2, 
data) ; 

xl = (xl + x2) * 0.5; 
yl = (yl + y2) * 0.5; 

z = aim_intensity (data->image, xl, yl) ; 
dzdx = aim_intensity (data->image, xl + 0.001, yl) - z; 
dzdy = aim_intensity (data->image, xl, yl + 0.001) - z; 
f_cartesian_to_polar (dzdx, dzdy) ; 
grad = fmag; 
if (grad == 0.0) 
grad = 0.0001; 
totalweight += length / grad; 
totalheight += z * (length / grad) ; 

} 

void measure_edge_contour difference (graphics data *data) 

{ 

snake s; 

int i , segments ; 
float *seg; 

double gap, pixsize, zoomfac; 

if (data->features->how_close >= 0.0) // if already computed 
return; 

boxes = data->f eatures->boxes; 
hi_x = data->f eatures->hi__x; 
hi_y = data->f eatures->hi_y; 




lo_x = data->features->lo_x; 
lo_y = data->f eatures->lo_y; 
box_si2e = data->f eatures->box_size; 
for (i = 0; contour_store [i] > -l.e6 || 

contour__store [i + 1] > -1.0e6; i += 2) { 

} 

segments = (i / 4) * 4; 

totallength = totalarea = totalsquares = 0.0; 
totalheight = totalweight = 0.0; 

sum_xdx = sum_ydy = sum_dx = sum_dy = sum_x = sum_y = 0; 
sum_wx = sum_wy = sum_x2 = sum_y2 = 0; 
pgf = 0; 

for (i = 0; i < segments; i += 4) { 
seg = ScContour_store [i] ; 
f ind_nearby_snakes (seg, data) ; 

} 

pixsize = (data->sem->pixel .X + data->sem->pixel .y) / 2; 

totalarea *= 0.5; 

totalsquares /= 3.0; 

if (totallength ==0.0) { 

printf ("Can' t match up contours with feature edges. \n"); 

data->f eatures->how_close = 999.0; 

return; 

} 

zoomfac = (double) data- >curpm- >zoom_l eve 1 / ZOOM_STEPS; 
gap = sqrt (totalsquares / totallength) * zoomfac; 
data->f eatures->how__close = gap * pixsize * 1000; 
gap = totalarea / totallength * zoomfac; 
printf ( 

"SEM/AIM difference: RMS = %.3f nM, Average = %.3f nM, Guess: 
%.4f \n", 

data->f eatures->how_close, gap * pixsize * 1000, 
totalheight / totalweight) ; 

} 

Subclaim 1-g -- excerpt from tracker. c -- Feb. 29, 2000 



Some of the code used in the mathematical algorithm that computes 
the next iteration's value for the above processing parameter is 
also used for other purposes and resides in tracker. c and plot.c. 

double aim_intensity (aimdata image, double X, double Y) 

{ 

int ny, ix, save_nx; 
double res; 

if (image == NULL) 

return (0 . 0) ; 
save_nx = nx; 
nx = image ->nx; 
ny = image ->ny; 

ix = (int) ( (Y - image->ll,y) / image->del .y) * nx + 

(int) ((X - image->ll.x) / image - >del .x) ; 
if (ix >= 0 && ix < nx * ny - nx - 1) 

res = interpolate (X, Y, image ->aimvalues + ix) ; 
else 



/ 



4 



res = 0.0; 
nx = save_nx; 
return (res) ; 



Subclaim 1-g excerpt from plot.c Feb. 29, 2000 



static compute_partial_derivatives (aimdata aim) 

{ 

Point_ptr ptr; 
int X, y; 
int anx, any; 
double sx, sy; 

anx = aim->nx; 

any = aim->ny; 

sx = 0.5 / aim->del.x; 

sy = 0.5 / aim->del.y; 

for (x = 0; x < anx; ++x) { 

ptr = aim->aimvalues + x; 

ptr->dzdy = 0; 

for (y = 1; y < any - 1; ++y) { 
ptr += anx; 

ptr->dzdy = (ptr[anx].z - ptr[-anx].z) * sy; 

} 

ptr += anx; 
ptr->dzdy =0; 

} 

for (y = 0; y < any; ++y) { 

ptr = aim->aimvalues + y * anx; 
ptr->dzdx = 0; 

for (x = 1; X < anx - 1; ++x) { 
++ ptr; 

ptr->dzdx = (ptr[l].z - ptr[-l] .z) * sx; 

} 

++ ptr; 

ptr->dzdx = 0; 

} 

} 

double interpolate (x, y, ptr) 
double X, y; 
Pointjptr ptr; 

{ 

Point_ptr yup; 

double delx, dely, fx, fy, yloz, yhiz, z; 
double sdlox, sdhix, sdx, sdloy, sdhiy, sdy; 
double qx, qy; 

yup = ptr + nx; 

/* Grid size */ 

delx = ptr[l] .X - ptr->x; /* uM */ 
dely = yup->y - ptr->y; 

/* Fraction of the way across the square */ 



fx = (x - ptr~>x) / delx; /* dimensionless */ 
fy = (y - ptr->y) / dely; 

/* Linear interpolation */ 

yloz = ptr->z * (1.0 - fx) + ptr[l].z * fx; /* Joules */ 

yhiz = yup->z * (1.0 - fx) + yup[l].z * fx; 

z = yloz * (1.0 - fy) + yhiz * fy; 

/* Quadratic correction */ 

sdlox = {ptr[l].dzdx - ptr[0].dzdx) / delx; /* J / uM^2 */ 
sdhix = (yup[l].dzdx - yup[0].dzdx) / delx; 
sdx = sdlox * (1.0 - fy) + sdhix * fy; /* J / uM^2 */ 

qx = sdx * { (fx - .5) * (fx - .5) - ,25) * delx * delx; 

sdloy = (yup[0].dzdy - ptr[0].dzdy) / dely; 
sdhiy = (yup[l].dzdy - ptr[l].dzdy) / dely; 
sdy = sdloy * (1.0 - fx) + sdhiy * fx; 

qy = sdy * ( (fy - .5) * (fy - .5) - .25) * dely * dely; 
z += qx + qy; 

/* Cubic correction - not yet implemented. */ 
return ( z ) ; 



